Phase structure of a spherical surface model 
on fixed connectivity meshes 



Hiroshi Koibuchi ^ 

^Department of Mechanical and Systems Engineering, Ibaraki National College of 
Technology, Nakane 866, Hitachinaka, Ibaraki 312-8508, Japan 



Abstract 

An elastic surface model is investigated by using the canonical Monte Carlo simula- 
tion technique on triangulated spherical meshes. The model undergoes a first-order 
collapsing transition and a continuous surface fluctuation transition. The shape of 
surfaces is maintained by a one-dimensional bending energy, which is defined on the 
mesh, and no two-dimensional bending energy is included in the Hamiltonian. 
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1 Introduction 



A considerable number of experimental studies have been conducted on the 
shape of membranes and their transformations jT]|2] . One of the interesting top- 
ics in biological physics is to understand the origins of such a variety of shapes 
statistical mechanically on the basis of membrane models [3P] . The curvature 
model of Helfrich, Polyakov, and Kleinert [^ll^lfT] is well known as a model for 
describing shapes of membranes and has long been used for understanding 
the fluctuation of membranes Pl9|10|ll|12|13j . From the two-dimensional dif- 
ferential geometrical view point, the curvature Hamiltonian is convenient to 
describe the mechanics of membranes 1131. 



However, the two-dimensional curvature Hamiltonian is not always necessary 
for providing the mechanical strength for the surface. In fact, it is well known 
that the cytoskeletal structures or the microtubules maintain shape of biolog- 
ical membranes 
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One-dimensional bending energy can serve as the Hamiltonian for a model of 
membranes. Skeleton models are defined by using the one-dimensional bending 
energy, which is defined on a sub-lattice of a triangulated lattice [15lll6fl7j . 
The compartmentalized structure constructed on the triangulated lattice is 
the sub-lattice and considered to be an origin of a variety of phases fTT] . 

The size of the sublattice is characterized by the total number n of vertices 
inside a compartment. As a consequence, the mechanical strength of the sur- 
face varies depending on n, because the compartment size is proportional to 
n and because the mechanical strength is given only by the sublattice. 

Therefore, it is interesting to see the dependence of the phase structure on 
n. The phase structure of the skeleton model is dependent not only on the 
bending rigidity h but also on the size n. The phase structure of compartmen- 
talized models at finite n was partly studied as mentioned above [T^16fl7j . 
On the other hand, the models are expected to be in the collapsed phase in 
the limit of n— >oo, because there is no source of the mechanical strength for 
the surface at n^oo. 

However, the phase structure in the limit of n — ^ is unknown and yet to 
be studied. The compartmentalized model in this limit is governed by one- 
dimensional bending energy and, for this reason the model is expected be 
different from the model with the standard two-dimensional bending energy 
defined on triangulated surfaces without the compartments. 

In this Letter, we study a triangulated surface model defined by Hamilto- 
nian that is a linear combination of the Gaussian bond potential and the 
one-dimensional bending energy. All the vertices are considered as junctions, 
which are considered to have a role for binding three one-dimensional chains. 
Consequently, the model is considered to be a compartmentalized model in 
the limit of n— >-0. Note also that the model in this Letter is allowed to self- 
intersect p^l9ll20] . and therefore the crumpled phase is expected to appear 
in the limit of 6— s-0, whereas the smooth phase is clearly expected in the limit 
of b^oo. 



2 Model and Monte Carlo technique 

Triangulated meshes are obtained from the icosahedron such that the bonds 
are divided into i pieces of the same length. Then we have the meshes of size 
N =10i'^ + 2, which include 12 vertices of coordination number g = 5 and, the 
remaining vertices are of coordination number q = 6. Figure [1] shows the mesh 
of size iV = 1442, which is given by £ = 12. The triangle surfaces are shown in 
the figure in order to visualize the mesh more clearly than that without the 
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triangle surfaces. 




Fig. 1. A triangulated mesh of size A^ = 1442, which is given by £ = 12. 

The Hamiltonian 5* is given by a linear combination of the Gaussian bond 
potential Si and the one- dimensional bending energy S2, which are defined by 

Si = Y.{X.-X,)\ 52 = E(l-t,-t,). (1) 



in 'S'l denotes the sum over bonds (ij), which connect the vertices i and 
j. In S2, tj is a unit tangential vector of the bond i. The symbol in 'S'2 is 
defined below. 

The pairing (ij) of the vectors tj and in S'2 is defined as follows: At the 
vertex of coordination number q = 6 such as O in Fig|2](a), we have three 
pairings 1—tAo " ^od, 1— t^o " ^oe, and 1—tco ■ ^of- The bending energy on 
the vertex of coordination number g = 5 such as O in Fig. MJo) is defined by 
five parings (l-t^io ■toc)/2, (l-tAo-toD)/2, (l-t^o •toD)/2, (l-tso ■to£;)/2, 
and (1— tco ■to£;)/2. Then, we effectively have 2.5 parings at the q = 5 vertices 
because of the factor 1/2; Z](ij) in >S'2 is defined by J2[ij) 1, where 1 = 1 at the 
vertices of g = 6 and 1 = 1/2 at the vertices of g = 5. Thus, we have 1 = , 
where Nb = 3N — 6 is the total number of bonds. 

The partition function of the model is defined by 

' N 

Z = I l[dX,exp[-SiX)], (2) 

1=1 

Six) = Si + bS2, 

where /' denotes that the center of the surface is fixed in the three-dimensional 
integrations /' Hili dXi, where Xi is the three dimensional position of the 
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Fig. 2. (a) A vertex of coordination number g = 6, and (b) a vertex of coordination 
number q = h. The bending energy in (a) is defined by three pairings 1— t^o ■ "toD, 
1— t^o ■ toE) and 1—tco ■ toF at the vertex O, and it is defined by the 2.5 pairings 
at the vertex O in (b) as stated in the text. 

vertex i. Because of the scale invariance of Z, Si is expected to be Si/N = 
3/2(A^-l)/Ar~3/2. 

The canonical Monte Carlo (MC) technique is used to simulate the multiple 
three-dimensional integrations in Z. The position X is shifted to X' = X+5X, 
where 5X is a position chosen randomly in a small sphere. The acceptance rate 
of the new position is about 50%. The radius of the small sphere is fixed at 
the beginning of the simulations. A random number sequence called Mersenne 
Twister pT] is used to the three-dimensional random shift and the Metropolis 
accept /reject in the MC simulations. 

The total number of Monte Carlo sweeps (MCS) at the region of the transition 
point after the thermalization MCS is 9 x 10*^ ~ 8 x 10^ for the A^ = 21162 and 
= 15212 surfaces, 6 x 10* ~ 5 x 10* for the = 10242 surface, 4 x 10* ~ 3 x 10* 
for the A^ = 7292 surface, and 3 x 10* ~ 2 x 10* for the A^ = 4842 and 
A^ = 2562 surfaces. Relatively small number of MCS is performed at non- 
transition region of h in each A^. 



3 Results 



Snapshots of surface of size A^ = 21162 are shown in FigsI31^a) and[3](b), which 
were obtained at 6 = 0.718 in the collapsed phase and in the smooth phase, 
respectively. Figures [3](c) and [3](d) show the surface sections in Figsj3]^a) and 
[3]^b). These four figures are shown in the same scale. The mean square size 
is about X2 = 79 in (a) and X2 = 145 in (b). 
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Fig. 3. Snapshot of surfaces of size N = 21162 obtained at 6 = 0.718 in (a) the 
coUapsed phase and (b) the smooth phase, (c) The section of the surface in (a), and 
(d) the section of the surface in (b). The mean square size is about = 79 m 
(a) and ^2 = 145 in (b). 

The mean square size is defined by 

X' = ^T.{X^^-X)\ X = Wx,, (3) 

i i 

where X is the center of mass of the surface. is expected to reflect the size 
or the shape of surfaces whenever the model has a smooth swollen phase and 
a collapsed phase. 

Figure m^a) shows versus b. The solid lines were obtained by the multihis- 
togram reweighting technique. The variation of X^ appears smooth against b, 
although it becomes rapid with increasing N. The variance Cx^ of X^ defined 
by 

Cx^ = ^{{x'-{X')y) (4) 

is plotted in Fig|l](b) against b. We clearly see in Cx^ an anomalous peak, 
which grows with increasing N. The anomalous peak seen in Cx^ represents 
a collapsing transition between the smooth swollen phase and the collapsed 
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Fig. 4. (a) The mean square size versus b, (b) the variance C^a of X'^ versus 
b, and (c) log-log plots of the peak values C^2^ against N. The solid lines in (a), 
(b) and the data C^f^ in (c) were obtained by the multihistogram reweighting 
technique. 

phase. 

In order to see the order of the transition, we plot the peak values C^2^ in 
FiglD^c) in a log-log scale against A^. The peak values 0^2^ and the statistical 
errors were obtained also by the multihistogram reweighting technique. The 
straight line in FigJl](c) was drawn by fitting the data to the scaling relation 

C^2^" oc N'', (5) 



where a is a scaling exponent. The fitting was done by using the data plotted 
in FigHl^c) excluding that of A^= 10242. Thus, we have 

a = 0.988 ±0.011, (6) 



which indicates that the collapsing transition is of first-order. The finite-size 
scaling (FSS) theory predicts that a transition is of first-order (second-order) 
if the exponent satisfies cr = 1 (a < 1). 

The Hausdorff dimension H is defined by ~ N'^/^ . We expect that ~ 2 
is satisfied in the smooth phase, whereas the value of H in the collapsed phase 
is unclear, because smoothly changes at the transition point as we see in 
FigHl^a). Therefore, in order to see the behavior of at the transition point 
more clearly, we plot the variation of X^ against MCS in Figsl5](a)-[5](i). The 
variations were obtained at 6 = 7.2, 6 = 7.22, 6 = 7.24 on the X = 10242 surface, 
at 6 = 7.16, 6 = 7.2, 6 = 7.22 on the X = 15212 surface, and at 6 = 7.16, 6 = 7.18, 
6 = 7.2 on the X = 21162 surface. 

We find from the figures that the value of X^ at the smooth phase is not so 
clearly separated from that of the collapsed phase at the transition point. In 
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fact, we can see a double peak structure only in the histogram h{X'^) of in 
FiglSj^h), although the double peaks are not so clear in the histogram, which 
is not depicted as a figure. No double peak structure was seen in h{X^) on the 
surfaces of iV< 15212. 

However, the mean value of at the smooth phase and that at the collapsed 
phase can be obtained from the series of X"^ in Figsl5](a)-[5]^i) by averaging 
between the lower bound X^;^ and the upper bound X^^^ assumed in each 
phase. Horizontal dashed lines in the figures denote X^^^^ and X^^^. 




Fig. 5. The variation of against MCS on the surface of size (a), (b), (c) A^ = 10242, 
(d), (e), (f) A^ = 15212, and (g), (h), (i) iV = 21162. The data were obtained at three 
distinct b close to the transition point in each A^, where those in (b), (e), and (h) 
are considered to be the ones obtained at the transition point. Horizontal dashed 
lines in the figures denote X^^^ and X^^x; which are shown in Tabled) 



The assumed values of X^™^ and X^^°Hn the collapsed phase and those of 
X^f ™° and X^|™° in the smooth phase are shown in Table [H The symbols col 
and smo denote the collapsed phase and the smooth phase, respectively. The 
values X^ in the collapsed phase on the surfaces of X = 21162, X = 15212, 
X = 10242, X = 7292 were respectively obtained at 6=7.16, 6 = 7.16, 6=7.2, 
6 = 7.16. On the other hand, those X^ in the smooth phase on the surfaces 
of X = 21162, X = 15212, X = 10242, X = 7292 were respectively obtained at 
6 = 7.2, 6 = 7.22, 6 = 7.24, 6 = 7.24. 

Figure [6] shows log- log plots of X^ versus X obtained in the smooth phase and 
in the collapsed phase. The straight lines were obtained by fitting the data to 
the relation X^ ~ X^/^, and we have the Hausdorff dimensions i^smo and i?coi 
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Table 1 

The assumed values of the lower bound X"^^ and the upper bound ^^ax 
obtaining the mean value in the collapsed phase close to the transition point, 
and those X^f™° and X^f™° in the smooth phase close to the transition point. 
6(col) and 6(smo) denote the bending rigidities where X'^ was obtained. 



N 


6(col) 


v2 col 
mill 


Y2 col 
^max 


6(smo) 


Y2 smo 
min 


Y2 smo 
^max 


21162 


7.16 


45 


100 


7.2 


115 


190 


15212 


7.16 


35 


85 


7.22 


95 


145 


10242 


7.2 


28 


65 


7.24 


73 
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7292 


7.16 
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41 


7.24 
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Fig. 6. Log-log plots of against N obtained in the smooth phase and in the 
collapsed phase close to the transition point. The straight lines were drawn by 
fitting the data to X^ ~ N'^^^ , where H is the Hausdorff dimension. 

respectively in the smooth phase and in the collapsed phase such that 

i^smo = 2.27 ± 0.29, H^oi = 2.29 ± 0.48. (7) 

The value of ifsmo is consistent to the expectation from the snapshot in 
Figsl3](b) and[3]^d). Moreover, we find from Hcoi in Eq.Q that the collapsed 
phase is considered to be physical, although Hcoi includes a large error. We 
must note that these values of H are dependent on the lower and the upper 
bounds X^^^ and -^^ax' therefore the results in Eq.([7]) are not so conclu- 
sive. Nevertheless, we feel that the phase transition of the model in this Letter 
is realistic. The physical condition Hco\ < 3 is expected to be obtained more 
conclusively by large scale simulations. 

Figure [Tl^a) shows the bending energy 82/ Nb versus b on the surface size 
N = 4842, = 10242, and = 21162. The reason for dividing ^2 by Nb is 
that in 5*2 of Eq.([T]) satisfies 1 = Nb as mentioned in the previous 
section. The slope of 82/ Nb becomes large with increasing N as expected. 
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Fig. 7. (a) The bending energy S2/NB versus 6, (b) the specific heat €$2 versus 
6, and (c) log-log plots of the peak values C™^^ against A^. The solid lines in (a), 
(b) and the data C^^^ in (c) were obtained by the multihistogram reweighting 
technique. The straight line in (c) was drawn by fitting the data to C^^^ oc . 

The specific heat defined by 
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S2 



52 



{S2)f) 



is plotted in Fig|7|(b). An anomalous peak can also be seen in Cs2 at the same 
transition point as that of the peak of C^a in FiglD^b). The peak values C^^^ 
are shown in Fig|7|(c) in a log-log scale against N . We draw in Figl7](c) the 
straight line which is obtained by the least squares fitting with the inverse 
statistical errors. The scaling relation is given by C^^^ oc A^'^, and we have v = 
0.325 ±0.019. Thus, we understand that the surface fiuctuation corresponding 
to the fiuctuation of 5*2 is a phase transition and is of second-order because of 
the argument of the FSS theory. 
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Fig. 8. (a) The bending energy S^/Nb versus b, (b) the specific heat C53 versus 
b, and (c) log-log plots of the peak values C^^^ against N. The solid lines in (a), 
(b) and the data C^^^^ in (c) were obtained by the multihistogram reweighting 
technique. The straight line in (c) was drawn by fitting the data to Cg^^^ oc A^'^. 
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The standard two-dimensional bending energy is defined by 5'3 = X](1— • rij), 
where iij is the unit normal vector of the triangle i. The bending energy 5*3 
is expected to refiect the surface fiuctuations, although it is not included in 
the Hamiltonian. Figure [Hl^a) shows S^/Nb versus 6, where the surface size 
is A^ = 4842, A^ = 10242, and Ar = 21162. The variance C53 = ^( {Sz-{S'i)f) 
defined by the expression similar to that of C^a in Eq.(jl]) is plotted in FigJHt^b), 
and the peaks C^^^ obtained by the the multihistogram reweighting technique 
are plotted against N in Fig|8]|^c) in a log-log scale. The straight line in Fig|8](c) 
was obtained by the least squares fitting, which was performed by using all the 
data in Fig|Hl(c). Thus, we have a scaling exponent ^ in the relation C^^^ccN^ 
such that /i = 0.307 ± 0.021. This result indicates that the surface fiuctuation 
transition is of second-order. 
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Fig. 9. The Gaussian bond potential Si/N versus b obtained on the toroidal surfaces 
of (a) A^ = 2562, iV = 4842, (b) A^ = 7292, iV = 10242, and (c) iV = 15212, iV = 21162. 
The expected relation S'i/A^~1.5 is satisfied. 



Finally, we plot in Figsl9](a)-[9](c) the Gaussian bond potential Si/N against 
b. As mentioned in the previous section, Si/N is expected to be S'i/A^~ 1.5 
because of the scale invariant property of the partition function and that of 5*2. 
This relation can always be used to check that the simulations were performed 
successfully. We see in the figures that the expected relation is satisfied. 



4 Summary and conclusions 



A triangulated surface model has been investigated by using the Monte Carlo 
simulation technique. Hamiltonian of the model is given by a linear combina- 
tion of the Gaussian bond potential and a one-dimensional bending energy. 
The model is considered to be obtained from a compartmentalized surface 
model in the limit of n — > 0, where n is the total number of vertices in a 
compartment and hence denotes the size of compartment. 

We have found that the model in this Letter undergoes a first-order coUaps- 
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ing transition and a second-order surface fluctuation transition. On the other 
hand, we know that the compartmentahzed model with the two-dimensional 
elasticity at the junctions undergoes a flrst-order surface fluctuation transition 
moreover a compartmentalized fluid surface model with the rigid junction 
also undergoes a flrst-order one |T7]. Therefore, we consider that the fluctu- 
ation of vertices inside the compartments strengthen the surface fluctuation 
transition in the n 7^ model. On the contrary, we have no vertices inside the 
compartments in the model of this Letter because of n = 0. The lack of vertex 
fluctuation is considered to soften the flrst-order surface fluctuation transition 
seen in the flnite n model. 



We should note that sufficiently small values of n implies that the compart- 
ment size is comparable to the bond length scale, which can arbitrarily be 
fixed due to the scale invariant property of the partition function. The size n 
is proportional to the area of a compartment, and hence the finite n implies 
that the corresponding compartment size is negligible compared to the surface 
size in the limit oi N ^ 00. The finite n also implies that the compartment 
size is sufficiently larger than the bond length scale. Thus, the model in this 
Letter is considered to be a compartmentalized model with sufficiently small 
compartment. 



The model in this Letter is allowed to self-intersect and hence phantom. A 
phantom surface model, which has a collapsing transition between the smooth 
phase and the collapsed phase, is considered to be realistic if the collapsed 
phase is physical. One of the criteria for such physical condition is given by 
H <3, where H is the Hausdorff dimension. Therefore, in order to see whether 
the condition is satisfied or not in our model, we obtained X"^ in the smooth 
phase and in the collapsed phase close to the transition point by averaging 
between X^^^ and X^^^ assumed in each phase. Thus, i^smo = 2.27(29) 
(smooth phase) and ifcoi = 2.29(48) (collapsed phase) were obtained, and then 
we found that the physical condition if < 3 is satisfied in the collapsed phase 
although Hcoi includes relatively large error. 



Meshwork models in has no vertex inside the compartments, which 

have finite size n. The phase structure of such meshwork model of finite n is 
considered to be dependent on the elasticity of junctions [23112^ . Therefore, it 
is interesting to study the dependence of the surface fiuctuation transition on 
n in the meshwork model, where the elasticity of junctions is identical to that 
in the model of this Letter. 
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